Short-time Dynamics of Percolation Observables 
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We consider the critical short-time evolution of magnetic and droplet-percolation order parameters 
for the Ising model in two and three dimensions, through Monte-Carlo simulations with the (local) 
heat-bath method. We find qualitatively different dynamic behaviors for the two types of order 
parameters. More precisely, we find that the percolation order parameter does not have a power-law 
behavior as encountered for the magnetization, but develops a scale (related to the relaxation time to 
equilibrium) in the Monte-Carlo time. We argue that this difference is due to the difficulty in forming 
large clusters at the early stages of the evolution. Our results show that, although the descriptions 
in terms of magnetic and percolation order parameters may be equivalent in the equilibrium regime, 
greater care must be taken to interprete percolation observables at short times. In particular, this 
concerns the attempts to describe the dynamics of the deconfinement phase transition in QCD using 
cluster observables. 

PACS numbers: ll.10.Cd, 11.15.Ha, 12.38.Aw, 14.70.Dj 



INTRODUCTION 

It is often useful to map a physical phase transition 
into the geometric problem of percolation, e.g. in order 
to gain a better understanding of how the transition is 
induced in the system (see e.g. 0). This mapping is ob- 
tained through a suitable definition of cluster in terms 
of the system's variables and parameters. In the case 
of the Ising model, such physical (droplet) clusters were 
introduced by Coniglio and Klein (based on the prescrip- 
tion by Kasteleyn and Fortuin) |2j and the mapping is 
well understood I3 ■ This holds also for the so-called n- 
vector models The correct cluster definition is of 
special interest in the description of more complex phase 
transitions, such as the quark-deconfinemcnt transition 
in finite-temperature QCD, for which the physical order 
parameter is not yet established. In this case there is still 
no satisfactory definition for a "droplet" cluster, even in 
the simpler pure-gauge case 

When equilibrium properties are investigated, the de- 
scriptions of the spin-model phase transition in terms 
of the magnetic order parameter or the percolation or- 
der parameter are indeed equivalent and one finds the 
same critical exponents (see, for example, and refer- 
ences therein). The same may or may not be true for 
the dynamical evolution of the different types of order 
parameters, although in principle one might expect to 
find equivalence for dynamic quantities as well, in par- 
ticular for the behavior at short times. We recall that 



the short-time behavior of the magnetic order parameter 
M is described by a scale-free expression, in terms of a 
power law (see for instance 0,11, El)- 

The study of the dynamic critical behavior of a per- 
colation order parameter with respect to the (Glauber) 
Monte-Carlo evolution might be of relevance for under- 
standing non-equilibrium effects in hot QCD, such as the 
effects due to heating and cooling of matter produced in 
heavy-ion collisions. Indeed, the possible connection be- 
tween the deconfinement transition in QCD and the per- 
colation phenomenon has received renewed attention 
in recent years || and the dynamics of cluster observables 
has been investigated using hysteresis methods . Note 
that the QCD phase transition is predicted to fall into the 
Ising universality class in the pure two-color [SU (2)] case 
and that the chiral phase transition in the two-flavor full- 
QCD case is expected to be in the universality class of 
the (continuous-spin) 4-vector model 0, 0, UJ. This 
motivates the connection between the percolation tran- 
sitions for spin models and the (dynamic) behavior at 
the QCD phase transition. We thus consider here the 
short-time dynamics of the two- and three-dimensional 
Ising model and focus on the dynamic critical behavior 
of the percolation order parameter. (Preliminary results 
of our study were presented in HE O.) Note that the 
dynamic behavior of cluster observables has been consid- 
ered for the droplet clusters in Ising and Potts models in 
various studies (see e.g. [13] and [lsj])- but mostly for the 
so-called cluster numbers n, — which denote the aver- 
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age number of clusters (per lattice site) containing s sites 
each — and not for the percolation order parameter. 

Here we investigate the dynamic critical behavior of 
the (zero-field) Ising model 



H 



-J ^ Si Sj 



(1) 



<ho> 



where J is positive and each Sj = ±1, for short Monte- 
Carlo times t, using the (local) heat-bath algorithm. We 
measure the magnetization M = (1/V) J2i$i an d we 
consider two different definitions for the percolation or- 
der parameter. Indeed, given a definition for a cluster on 
the lattice, one may consider as the order parameter in 
percolation theory the stress of the percolating cluster, 
defined for each configuration by the relative volume of 
the infinite cluster. (In particular, in the droplet picture 
for the Ising model, this would give zero for temperatures 
above the critical one.) On a finite lattice, one might con- 
sider the volume of the infinite cluster to be the one of 
the spanning cluster, or zero when there is no percolating 
cluster (i.e. above the critical temperature). Note that a 
spanning cluster is a set of spins connected from the first 
to the last row of the lattice in at least one of its space di- 
rections. Alternatively, one can also consider the relative 
volume of the largest cluster (see e.g. Chapter 3 in £|). 
We denote these two definitions for the percolation order 
parameter respectively by Q (in the case of the spanning 
cluster) and CI' (in the case of the largest cluster) . 

In this work we compare the behaviors of the two types 
of quantities just described (i.e., percolation or magnetic 
order parameters) as functions of t. As shown below, we 
find that whereas the magnetic order parameter displays 
a power-law increase with t, the data for the percolation 
order parameters CI and CI' are well fitted by a diffusion 
and by a growth and nucleation process respectively. In 
both cases the time scales — respectively called t and 
t' — are related to the relaxation time to equilibrium. 
Indeed, we show that r and r' diverge as L z when the 
lattice side L tends to infinity, where z « 2 is the dynamic 
critical exponent of the heat-bath algorithm 



SHORT-TIME (MONTE CARLO) DYNAMICS 

Using renormalization-group theory, it can be shown 
0, that the early time evolution of an order parame- 
ter (e.g. the magnetization M) already displays universal 
critical behavior, given by 



M(t,e,m ) = b- p/ " 'M{tb' z ,eb 1 /" \m b x °) 



(2) 



where tuq is the initial magnetization, e = (T — T c )/T c , 
Ai is a universal function and b is a scale factor, which 
can be taken equal to t x l z . We thus expect for T = T c 
and small mo a power-law behavior at early times 



with 8 = (xq — [3/v)/z . In principle, we would assume 
that the two percolation order parameters CI or CI' defined 
above should have a similar behavior. 

The heat-bath dynamics consists in choosing the two 
possible directions of each Ising spin according to the ex- 
act conditional probability given by its nearest neighbors. 
Each spin Sj is chosen "up" or "down" respectively with 
probability pi or 1 — , where 



I 



Pi = 



l + exp(-2/3E,^ 



(4) 



and the sum is over nearest neighbors of Si. After a 
certain number of iterations the spin configuration obeys 
the Boltzmann distribution. In the heat-bath method, 
since the updates are local, this transient time becomes 
considerably large at criticality. 

The Fortuin-Kasteleyn clusters are obtained from the 
Ising-model Hamiltonian by writing the partition func- 
tion as 



(5) 



m m i (id) 



M(t) £ _ ~ mot 6 , 



(3) 



where pij = 1 — exp(— 2 J/3) is the probability of having 
a link between two nearest-neighbor sites of equal spin 
value. This link is represented by mj and determines the 
clusters that will be associated with percolation at the 
critical temperature [2l|. Note that the above defined 
clusters are used in the Swendsen-Wang algorithm 
to perform global moves in which the spins in a cluster 
are flipped together. Here we only use these clusters to 
calculate percolation observables, whereas the dynamics 
is given by local heat-bath updates, as described above. 



NUMERICAL RESULTS 

In order to study the short-time dynamics we simu- 
late at T = T c and force the system to have an initial 
magnetization mo. We let the system evolve in time and 
look for power-law behavior of the order parameters M, 
CI and Cl' as functions of the (Monte-Carlo) time. Each 
temporal sequence is generated from a different random 
seed, i.e. each sequence has a different initial spin con- 
figuration. The time history is then obtained from an 
average over all the generated sequences. 

We have studied the two- and three-dimensional cases, 
performing Monte Carlo simulations respectively with 
50,000 and 40,000 seeds, for several initial magnetizations 
mo and several lattice volumes, using the heat-bath algo- 
rithm. We consider periodic boundary conditions. Note 
that in order to compare the percolation order parame- 
ters Cl and Cl' to M we consider the volumes of clusters of 
"up" spins with positive sign and volumes of clusters of 
"down" spins with negative sign when taking the average 
over the seeds. 
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We obtain that a power-law fit works very well for M, 
yielding the literature value [2^| for the exponent 6. How- 
ever, as can be seen in Figs.nand|21 the percolation order 
parameters and O' do not show a power-law behavior, 
being consistent with exponential behaviors in terms of 
t/r, thus having r as a time scale. As verified below, 
the exponential behaviors are different for Q and f2', but 
t is in both cases directly related to the relaxation time 
to equilibrium. We thus find the surprising result that 
although the equilibrium behaviors of M, £1 and Q' are 
equivalent, the different types of order parameters show 
qualitatively different dynamic critical behavior. 
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TABLE I: Fits of fi(t) to the form A exp(-r/t). Data for the 
case mo = 0.02. The fit intervals are chosen as shown in the 
second column, whereas the errors in the fit parameters are 
adjusted to account for the fluctuations in these parameters 
for slightly different fit intervals. 



Volume 


At 


A 


T 


2/7 r 

X /d.o.f. 


100 2 


18-83 


0.061(1) 


85.5(9) 


0.97 


125 2 


28-110 


0.062(1) 


134(1) 


0.98 


150 2 


30-150 


0.063(1) 


192(2) 


0.91 


200 2 


41-240 


0.061(1) 


334(3) 


0.69 


24 3 


3-100 


0.0366(3) 


5.9(3) 


0.57 


32 3 


10-120 


0.0397(5) 


11.7(6) 


0.42 


48 3 


15-120 


0.0446(4) 


28.4(4) 


0.58 


64 3 


20-130 


0.0467(3) 


51.6(3) 


0.68 



and three-dimensional cases by 

fi(t) = A exp(-r/t) . 



(6) 



as can be seen from Table [I] On the other hand, the 
behavior of fi' is given by 



n'(t) = B{1 - exp [-(*/•/)"]} 
as shown in Table CU 



(7) 



FIG. 1: Plot of the early time evolution of the magnetic 
(M) and percolation (fi, fi') order parameters for the two- 
dimensional case. Data are shown for mo = 0.02 and L = 200. 
Note the logarithmic scale on both axes. 
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FIG. 2: Plot of the early time evolution of the magnetic 
(M) and percolation (fi, Q') order parameters for the three- 
dimensional case. Data are shown for mo = 0.02 and L — 64. 
Note the logarithmic scale on both axes. 



More precisely, the behavior of fl is given in the two- 



TABLE II: Fits of fi'(t) to the form 5(1 — exp[-(t/-r' )"])• 
Data for the case mo = 0.02. The fit intervals are chosen 
as shown in the second column, whereas the errors in the 
fit parameters are adjusted to account for the fluctuations in 
these parameters for slightly different fit intervals. 



Volume 


At 


B 


r 


V 


X 2 /d.o.f. 


100 2 


7-100 


0.104(5) 


82(1) 


0.71(1) 


0.35 


125 2 


8-100 


0.126(5) 


168(14) 


0.70(1) 


0.23 


150 2 


7-100 


0.115(4) 


179(13) 


0.72(1) 


0.17 


200 2 


10-100 


0.20(1) 


380(32) 


0.73(1) 


0.17 


24 3 


3-100 


0.0349(1) 


8.0(1) 


0.82(1) 


0.42 


32 3 


5-120 


0.0372(2) 


15.4(3) 


0.83(2) 


0.50 


48 3 


10-130 


0.0394(3) 


34.4(4) 


0.90(1) 


0.48 


64 3 


15-150 


0.0402(5) 


60.5(4) 


0.93(1) 


0.56 



Fits of the data to the above forms are shown (for the 
largest lattices considered) in Figs. and ^respectively 
for the 2d and 3d cases. We see that the departure of 
the percolation order parameters from the power-law be- 
havior of the magnetization remains even long after the 
so-called microscopic time 0. In other words, neither Q, 
nor fr!' can be fitted to a power law at short times. In 
fact, the forms proposed for and W in Eqs. (© and @ 
are well fitted by the data in two and three dimensions, 
for all volumes considered. (We note however that the 
two-dimensional data for Q were fitted in [lr3 | considering 
also a second exponential term, adding a third parameter 
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to this fit.) The above data for r and r' can be fitted 
to the form L 2 , in the two- and three-dimensional cases 
separately. In all four cases we find z « 2. The two 
time scales can therefore be associated with the relax- 
ation time to equilibrium. 
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FIG. 3: Data for the 2d case and fits of Q. (lowest curve) and 
Q! (middle curve) according to the forms in Eqs. @ and Q 
respectively. The magnetization (top curve) is also shown, 
for comparison. Data and fits are shown for mo = 0.02 and 
L = 200. Error bars correspond to one standard deviation. 
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FIG. 4: Data for the 3d case and fits of SI (lowest curve) and 
Si' (middle curve) according to the forms in Eqs. © and Q 
respectively. The magnetization (top curve) is also shown, 
for comparison. Data and fits are shown for mo = 0.02 and 
L = 64. Error bars correspond to one standard deviation. 

Note that the behavior of SI (in two and three di- 
mensions) is a factor \ft times the solution of a one- 
dimensional diffusion equation at a fixed position in space 
[24]|. In fact, the expression for the concentration of diffu- 
sion material due to a change in chemical potential along 
a direction 1 at a position xo and time t is given by 
(C/y/t) exp(— x^/ADt) , where D is the diffusion coeffi- 
cient. Comparing this to the form for S7 in Eq. JJJJ, we 



see that if we multiply the percolation order parameter 
in this case by \[t we get the diffusive form, identifying 
Xq/AD with r. Of course, the factor y/i corresponds to 
the diffusion length for a random walk at time t. More- 
over, the fixed position xq is proportional to y/r, i.e. 
roughly the length of the lattice (since r as L z as men- 
tioned above) . This connects the occurrence and strength 
of percolation at time t to the probability that a random 
walk will reach a length L after t steps. 

The behavior of Q! corresponds to the volume fraction 
of particles in a process of nucleation and growth [2^ . 
This behavior is observed in the dynamics of weak first- 
order phase transitions as seen e.g. in Ref. [25J. In that 
reference, the critical dynamics of a scalar field quenched 
to a metastable state was investigated with a model-A 
Langevin equation. In this context, the result for SI' is 
particularly illuminating in the sense that for the same 
underlying dynamics (model A, in this case of the heath- 
bath dynamics), the short-time behavior clearly depends 
on the specific observable and its initial configuration. It 
would be also interesting to check if the exponent 77 goes 
to 1 in the limit of very large L in three dimensions (as 
seems to be suggested by the data in Tabled or to relate 
the exponents in two and three dimensions to the fractal 
dimensions associated with the percolating clusters in the 
two cases. 

Regarding the so-called percolation cumulant or per- 
colation probability — taken as 1 if there is percolation 
and if there is not — we obtain that this quantity does 
not show a power-law behavior in time, contrary to what 
is observed for the Binder cumulant . We find that the 
percolation cumulant is also described by an exponential 
exp(— r/t), with the prefactor \ft in two dimensions and 
£-0.34 j n three dimensions. 



CONCLUSIONS 

We have investigated numerically the critical heat-bath 
dynamics for magnetic and percolation order parameters 
in the Ising model at short Monte Carlo times, starting 
from a small magnetization mo- From our results we see 
that although the equilibrium behaviors of the magneti- 
zation M and of the percolation order parameters SI and 
£7' are equivalent, the two types of order parameters show 
qualitatively different dynamic critical behavior at short 
times. Indeed, whereas the magnetic order parameter 
M shows a power-law behavior with the exponent 0, one 
finds that SI and SI' have a time scale, given respectively 
by r in Eq. © and r' in Eq. Q). This time scale seems 
to be related to the relaxation time of the algorithm used 
for thermalization. 

The short-time behaviors of SI and SI' are well de- 
scribed respectively by diffusion and by growth and nu- 
cleation processes, probably related to the difficulty in 
forming a percolating cluster at the early stages of the 
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simulation. 
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